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ABSTRACT 



The magnetorotational instability (MRI) plays a crucial role for cosmic struc- 
ture formation by enabling turbulence in Keplerian disks which would be other- 
wise hydrodynamically stable. With particular focus on MRI experiments with 
liquid metals, which have small magnetic Prandtl numbers, it has been shown 
that the helical version of this instability (HMRI) has a scaling behaviour that is 
quite different from that of the standard MRI (SMRI). We discuss the relation 
of HMRI to SMRI by exploring various parameter dependencies. We identify 
the mechanism of transfer of instability between modes through a spectral ex- 
ceptional point that explains both the transition from a stationary instability 
(SMRI) to an unstable travelling wave (HMRI) and the excitation of HMRI in 
the inductionless limit. For certain parameter regions we find new islands of the 
HMRI. 

Subject headings: accretion, accretion disks; instabilities; MHD; turbulence 
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Introduction 



The magnetorotational instability (MRI) (IBalbus 20091 ) is considered as the main 
candidate to solve the long-standing puzzle of how stars and black holes are fed by the 
accretion disks surrounding them. The central problem is that these accretion disks 
typically rotate according to Kepler's law, f2(r) ~ r ~ 3 / 2 5 w hich results in an angular 
momentum r 2 f2(r) ~ r 1 / 2 . Hence, they fulfill Rayleigh's criterion stating that rotating 
flows with radially increasing angular momentum are hydrodynamically stable, at least in 
the linear sense. Such stable, non-turbulent disks would not allow the outward directed 
angular momentum transport that is necessary for the infalling disk matter to accrete into 
the central object. 



In their seminal paper of 1991 (IBalbus fc Hawley 19911 ) Balbus and Hawley had 



highlighted the key role of the magnetorotational instability (MRI) in explaining turbulence 

and angular momentum transport in accretion disks around stars and black holes. They 

had shown that a weak, externally applied magnetic field serves only as a trigger for the 

instability that actually taps into the rotational energy of t he flow. This is quite in contrast 

I 

to current-induced instabilities, e.g. the Tayler instability (ITayler 19731 ). which draw their 
energy (at least partly) from the electric currents in the fluid. 

Soon after the paper by Balbus and Hawley it became clear that the principle 
mechanism of t he MRI had already been revealed three decades earlier by Velikhov 



(jVelikhov 19591 ) and Chandrasekhar (IChandrasekhar 19601 ) . Actually, they had investigated 



the destabilizing action of an external magnetic field for the classical Taylor-Couette (TC) 
flow between two concentric, rotating cylinders rather than for Keplerian rotation profile. 
This is, however, not a crucial difference since a TC flow can be made very close to a 
Keplerian one simply by adjusting the ratio of rotation rates of the inner and the outer 
cylinder. 
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The MRI in flows between rotating walls has attracted renewed interest during the last 
decade, mainly mot i vated by the increas ing efforts to investigate MRI in the laboratory 
f lRosner et al. 20041 ; 



Stefani et al. 2008al ). A first inter esting experiment al result was 



obtained in a spherical Couette flow of liquid sodium (ISisan et al. 20041 ). The authors 
observed correlated modes of velocity and magnetic field perturbation in a parameter region 
which is quite typical for MRI. It must be noted, however, that the background state in this 
spherical Couette experiments was already fully turbulent, so that the original goal that the 
MRI would destabilize an otherwise stable flow was not met. At Princeton University work 
is going on to identify MRI in a TC experiment with liquid gallium, and first encouraging 
results, in cluding the observation of non-axisym metric Magneto- Coriolis waves, have been 



obtained (INornberg 2008 



Nornberg et al. 20091 ). 



Both experiments had been designed to investigate the standard version of MRI (SMRI) 
with only a vertical magnetic field being applied. In this case, the azimuthal magnetic field 
(which is an essential ingredient of the MRI mode) must be produced from the vertical field 
by induction effects, which are proportional to the magnetic Reynolds number (Rm) of the 
flow. Rm, in turn, is proportional to the hydrodynamic Reynolds number according to 
Rm = PmRe, where the magnetic Prandtl number Pm = vjr\ is the ratio of viscosity v to 
magnetic diffusivity rj = \j\i^a. For liquid metals Pm is typically in the range 10~ 6 ...10~ 5 . 
Therefore, in order to achieve Rm ~ 1, we need Re ~ 10 5 ...10 6 , and wall-constrained flows 
(in contrast to wall-free Keplerian flows) with such hi gh Re are usua lly turbulent, whatever 



the linear stability analysis might tell (see, however, ( IJi et al. 20061 ) ). This is the point 



which makes SMRI experiments, and their interpretation, so cumbersome. 

One might ask, however, why not to substitute the induction of the necessary azimuthal 
magnetic field component of the MRI mode by simply exte rnally applying this compone nt 



as a part of the base configuration. Indeed, it was shown (IHollerbach fc Riidiger 2005 
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Rudiger et al. 20051 ) that the resulting "helical MRP (HMRI), as we now call it, is then 
possible at far smaller Reynolds numbers and magnetic field amplitudes than SMRI, making 
HMRI an ideal playground for liquid metal experiments. 

First experimental evidence for HMRI was obtained in 2006 at the liquid metal facility 
PROMISE (Potsdam POssendorf Magnetic /notability Experiment) which is basically a 



Taylor- Couette (TC) cell made of concentric rotati ng copper walls, fi 



eutectic which is liquid at room tempera tures). In (jStefani et al. 2006 



l ed with GalnSn (a 



Stefani et al. 2007 



Riidiger et al. 2006 



Stefani et al. 2008bl ) it was shown that the HMRI travelling wave 
appears only in the predicted finite window of the magnetic field intensity, with a frequency 
of the travelling wave that was also in good accordance with numerical simulations. Results 
of a significantly improved experiment (P ROMISE 2) with strongly reduced Ekman pumping 



at the end-caps were published recently ( jStefani et al. 2009a 



Stefani et al. 2009bh 



The connection of SMR! 



Riidiger fc Hollerbach 2 007 



and HMRI is pre s ently under intense deba t e (|Liu et al. 2006 ; 



Priede et al. 2007 



Lakhin fc Velikhov 2007: 



Liu et al. 2007 



Szklarski 2007 



Riidiger & Schultz 2008 



Liu 20091: 



Priede fc Gerbeth 2009). The first 



essential point to note here is that HMRI and SMRI are connected. Indeed, Fig. 1 in 



( iHollerbach fc Rudiger 20051 ) shows that there is a continuous and monotonic transition 



from HMRI to SMRI when Re and the magnetic field strength are increased simultaneously. 
A second remarkable property of HMRI for small Pm (which has been coined 



"inductionless MRP), was clearly worked out in (IPriede et al. 20071 ). It is the apparent 
paradox that a magnetic field is able to trigger an instability although the total energy 
dissipation of the system is larger than without this field. 

The relevan ce of HMRI for Keplerian flows has been seriously put into question in 



( ILiu et al. 20061 ). Using a local WKB analysis in the small-gap approximation, the authors 



had shown that the HMRI works only for comparably steep rotation profiles (i.e. slightly 
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above the Rayleigh line) and disappears for p rofiles as flat as the Keple rian one. This result 



has bee n confirmed by Lakhin an d Velikhov (jLakhin fc Velikhov 20071 ) and Riidiger and 



Schultz flRiidiger fc Schultz 20081 ). 



However, this disappointing result was soon relativized in (jRudiger fc Hollerbach 20071 ) 
by solving the global eigenvalue equation for HMRI with electrical boundary conditions. It 
turned out that HMRI re-appears again for Keplerian flows provided that at least one radial 
boundary is highly conducting. A similar discrepancy between local and global resu lts is 



well known for the so-called stratorotational instability (SRI) (IDubrulle et al. 20051 ) for 



which the existe nce of reflecting boundaries appears necessary for the instability to work 



(jUmurhan 20061 ). This artificial demand is of course a much stronger argument against the 
working of SRI than the necessity of one conducting boundary is for the working of HMRI: 
considering, i.e., the colder outer parts of accretion dis ks, then the inner part can indeed be 



considered as a good conductor (IBalbus fc Henri 20081 ). 



Another argument that has been put forward against the relevance of HMRI for thin 
accretion disks is the necessity for a large ratio of toroidal to poloidal magnetic fields 



flLiu 20081 ). 



A further complication for applying HMRI to the real world is the fact that it appears 
in form of a travelling wave. The crucial point here is that monochromatic waves are 
typically not able to fulfill the axial boundary conditions at the ends of the considered 
region. To fulfill them, one has to consider wave packets. Only wave packets with vanishing 
group velocity will remain in the finite length system. Typically, the onset of this absolute 
instability, characterized by a zero growth rate and a zero group velocity, is harder to 
achieve than the convective instability of a monochromatic wave with zero growth rate. A 
comprehensi ve analysis of the relatio n of convective and absolute instability for HMRI can 



be found in (IPriede &: Gerbeth 20091 ). From the extrapolation of the results of this paper it 
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seems that Keplerian rotation profiles (with conducting boundaries) are indeed absolutely 
HMRI-unstable, but a final solution to this puzzle is still elusive. 

In the present paper, we step back from those important consideration of absolute 
and global instabilities and focus again on the local WKB method b y considering th e 



dispersion relation of MRI which had been deriv ed and analyzed in (ILiu et al. 2006 



Lakhin fc Velikhov 2007 



Riidiger fe Schultz 20081 ). In spite of these former investigations, 



we feel that some points still need further clarification. This concerns a careful application of 
the Bilharz stability criterion as well as some further parameter dependencies, in particular 
the dependence for small but finite magnetic Prandtl numbers. It also concerns the question 
in which sense the HMRI can be consid ered as a dissipation induced instability which is 



quite common in many areas of physics (jKrechetnikov fc Marsden 2007 : 



Kirillov 20070 . 



To make the paper self-contained, we will start with a re-derivation of the dispersion 
relation in two forms which explicitly contain the relevant frequencies or the dimensionless 
parameters, respectively. 

Then we will study the peculiar relation of SMRI and HMRI. As a main result of this 
paper we will describe in detail the mechanism of transition from SMRI to HMRI through 
a spectral exceptional point which appears at finite but small Pm. This provides a natural 
explanation for the continuous and monotonic connection between SMRI (a destabilized 
slow magneto-Coriolis wave) and HMRI (a weakly destabilized inertial oscillation). In 
addition to this, for high Reynolds numbers we will identify a second scenario for HMRI 
which leads to new islands of instability at small but finite values of Pm. 
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2. Mathematical setting 

In order to make the paper self-contained, we will re-derive in this section the 
dispersion relation for HMRI, including viscosity and resistivity effects. Note that this 
dispersion relation was derived in various f orms and approxima tions by a number of authors 



flLiu et al. 2007 



Lakhin k Velikhov 2007 



Riidiger et al. 20081 ). 



Th e standard set of non-linear equation s of dissipative incompressible magnetohydrody 



namics ( 


Ji et al. 2001 


■ 


Goodman k Ji 2002: 


Noeuchi et al. 2002 


Rudieer k Schultz 2008 


>) consists of the Navier-Stokes equation : 



Lakhin k Velikhov 2007 



<9u 

~dt 



-vu + — 

P \ 2/xq 



1 



Hop 



[u- V)u = 

and of the induction equation for the magnetic field B 

Vx(uxB)+j]V 2 B, 



B • V)B + z/V 2 u, 



dB 

~dt 



(2) 



where p is the pressure, p = const the density, u = const the kinematic viscosity, 
n = (/iocr) -1 the magnetic diffusivity, a the conductivity of the fluid, and po the magnetic 
permeability of free space. Additionally, the mass continuity equation for incompressible 
flows and the solenoidal condition for the magnetic induction yield 



V • u = 0, V ■ B = 0. 



(3) 



We consider the rotational fluid flow in the gap between the radii R\ and R2 > Ri, 
with an imposed magnetic field sustained by currents external to the fluid. The latter is 
important in order to distinguish the MRI from other instabilities (i.e. the Tayler instability 
for which electric currents are applied to the fluid). Introducing the cylindrical coordinates 
(R, (j), z) we consider the stability of a steady-state background liquid flow with the angular 
velocity profile Q(R) in helical background magnetic field (a magnetized Taylor-Couette 
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flow) 

u = RQ(R) e^, p = p Q (R), B = B°(i2)e + B° z e z 
with the azimuthal component 



B°AR) 



2ixR' 



(4) 



(5) 



which can be thought as being produced by an axial current I. The angular velocity profile 
of the background TC flow is 



n(R) = a + 



R r 



(6) 



where a and h are arbitrary constants as in Taylor- Couette experiments (IWendl 19991 ) . The 
centrifugal acc eleration of the background flow (j5D is compensated by the pressure gradient 
flJi et al. 20011 ) 



RQ 



, = I dpp 
pdR' 



(7) 



2.1. Linearization with respect to axisymmetric perturbations 

Throughout the paper we will restrict our interest to axisymmetric perturbations 
u' = u'(R, z), B' = B'(i2, z), and p' = p'(R, z) about the stationary solution (jl])-([7j), keeping 
in mind that for strongly dominant azimuthal m agnetic fields also non-axisymmetric 



perturbations are possible (IHollerbach et al. 2009i ). 
With the notation 

D^dnd^ + dl D 2 = did R + dl (8 

where the differential operators are defined in flAip . the general linearized equations flA2[) 
derived in the 
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Appendix are simplified in the assumption of axisymmetric perturbations to 



(d t - vD 1 )v! R - 2Qu' 



(d t - vDxju'j, + —u' R 

(d t - vD 2 )u' z 
(dt-rjDjB'x 
(ft-i^iJBj 
(d t -rjD 2 )B' z 
d 7 u' 



d z B' z 



d z B' 



1 

P 
Bl 

PoP 
1 

P 



B%u' R , 



d R p' + - (B° z d R B' z + Bld R B'^) 
Po 



f 



d zP ' + - (B%B' Z + B\d z B^) 
Po 



9R° 



B%u' + — V fl + (Rd R Q)B' R , 



R 



-B°aU 



-9rU r , 
-d R B' R , 



+ — I B%B' R - 
Hop V R 



3 

PoP 



■d z B' z , 



(9) 



where the squared epicyclic frequency k is defined as 



Following the approach of (IGoodman fc Ji 2002 



(10) 



Liu et al. 20061 ) we act on the first of 



equations (J9J) by the operator d R and on the third one by the operator d z . Summing the 
results, taking into account that 



d R D x = D 2 d R , 9^ = 0, d R (B°B') 



m 

R 



+ B° 6 d R B' 



(11) 



and using the last two equations of ([9]) yields 



2Qd R u 4 



1 

P 



-Do 



1 



V + ~B° Z B 
Po 



1 1 

PPo 



{dR Bld R + dlBl)B',-^ p d R {^B', 



- l -D 2 
P 



p> + -(B° Z B' Z + B\B\ 
Po 



Therefore, we extend the identity obtained in (IGoodman &: Ji 2002l ) to the case I?9 ^ 



2 



PoP 



R 



(12) 



Po 



(13) 
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On the other hand, using (fTTl) we transform the first of the equations (Q into 



(d t -pD 1 )u' R -2Qu , (j> = -d. 



R- 



P 



p' + - (B°B Z + tiffi) 
Po 



+ 



2B<: 



-B' 



(14) 



p p V " R 

Acting on both sides of the equation (j!4p by the operator D\ and taking into account the 
identity ffT3l) and 



D.dn = (d R d R + d 2 z )d R = d R (d R d R + d 2 z ) = d R D- 



(15) 



we get 



D x {d t -vD x )u' R - 2QD lU ' 



d R D 2 - 



+ D 1 —[B%B'-'^B' 



1 

Hop 



p' + - (B°B' Z + BlB',) 
PL Po 

2d R d R ( nu' 6 - — ^-B' 6 )+D 1 —( B<!d,B'„ - -J-Bl 



R 



Pop R v pop V R 
Rearranging the terms and using the definition of the operator D\ yields 



(16) 



D 1 (d t -uD 1 )u' R -2Qd 2 z u f (j) = 2d R d R —^-B'+—B z D l d z B' R -—D 1 ^B i 



1 2Bl 

1 n <t> r>i 



Pop R 



PoP 



Pop R 



1 i 9/?° 

^B° z D 1 d z B' R -—^td*Bl 



pop ~ " pop R 

Therefore, we have separated the equations for u' R , u', and B' R , B', from the others in 



(17) 



(d t - vD x )D x v! R - 2VLd 2 z u' 



(d t - uDx)^ + — u' R 



1 i 2B° 

l -BlD ldz B' R -^^dlB' 



PoP 
PoP 



Pop R 



d z B' 



(dt-vD^B^ 



B%u' R , 



25 



B °' dzU * + R U ' R + ( RdRn)B *- 



Note that after introducing the stream functions for the poloidal components 
u' R = d z ip, u' z = -d R ip, B' R = d z ip, B' z = -d R ip, 



(19) 



-12- 



the equations (fTBj) extend the inviscid equations of (ILiu et al. 20061 ) to the case v ^ 0. 



We can rewrite (ITS]) in the form of the operator matrix equation d t E£,' = H£', where 

/ 



E 



^ D 1 ^ 



10 
10 
1 



V 



H 



J 



y z flop 1 z 





B"i.h T)D X 
B% Rd R Q 



2B° 



2B l B 2 \ 
HopR * 

POP 





(20) 



/ 



The resulting multiparameter family of operator matrices equipped with boundary 
conditions can be investigated by numerical or perturbative ( Kirillov 201 ol ) methods. In the 
following we use the local WKB approximation. 



2.2. Local WKB approximation 



We choose a fiducia l point (Rq, Zq), around which we perform the local stability analysis 



(IPessah fc Psaltis 20051 ). We expand all the background quantities in Taylor series around 



(Rq, zq) and retain only the zeroth order in terms of the local coordinates R = R — Rq and 
z = z — z to obtain the operator matrix equation with the constant coefficients 



(21) 



with 



-En 



' D\ ^ 

10 

10 

1 



Hq 



v(Dl 



0\2 

















POP 1 





2B' 



popRo 

Pop 




(22) 



13 



where 



dR 



R=Ro , 



Ro Rq 
(23) 



Equation fl2Tl) is a linear PDE with the constant coefficients in the local variables (R, z) 
for the perturbed quantities This is a good approximation as long as the variations R 
and z are small in comp arison with the charact eristic length scales in the radial and vertical 



directions, respectively (IPessah fc Psaltis 20051 ). A solution to the equation (I2TI) has the 
form of a plane wave 



£' = £exp (jt + ik R R + ik z z), | = {u R ,u<j>, Br, , 

where £ is a vector of constant coefficients. 

Introducing the total wave number k 2 = k 2 + fc^ and denoting a = k z /k, we find 

ik R 1 



(24) 



-D° exp (jt + ik R R + ik z z) 



k H — — — — — 2 ) exp [p/t + ik R R + ik z z). (25) 



i? -^o, 

In the WKB approximation we restrict the analysis to the modes with the wave numbers 
satisfying k R R ^> 1 which allows us to neglect the terms ^ — in (1251) . In view of this, 
after substitution of (1241) into equation fT2Tl) . we arrive at the matrix eigenvalue problem 



(#- 7 £)f=0, 



(26) 



with £ as a unit matrix and H = —dia,g(uj u ,uj u ,uj v ,uj v ) + Hi + H2, where u u = vk 2 and 
uj v = r/k 2 are the viscous and resistive frequencies, 



Hi 









1 














1 


PoP 














PoP 









(27) 
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/ 



Ho 



V 

and the Alfven frequencies are 






2Q a 2 





- 2w ^ 


oo r? dSl 1 
ZU -K ^\ r=Rq 




























73 dQ 1 
^0 dR\R=R 






(28) 



0\2 



0\2 



PoP 



PopR 



2 ■ 



(29) 



Note that the matrix — diagfa;,,, uj u , u v , ujrj) + Hi has two double eigenvalues related to 
the damped Alfven modes ( Nornberg et al. 2OO9I ) 



7l,2 



u v + uj v 




u„ — uo„ 



— UJ 



A- 



73,4 



U v - Ur 



u . 



(30) 



When ua s = 



dn I 

' dR\R=R 



0, the eige nvalues of the matrix Hi + H2 correspond to the 



Alfven-inertial or Magneto-Coriolis waves (ILehnert 1954 ) 



71,2 = n/^A 



0, 13,4 = A 



(31) 



Figure [T] demonstrates how rotation le ads to the spli tting of plane Alfven waves into the 



fast and slow Magneto-Coriolis waves (ILehnert 19541 ) . The system with purely imaginary 



eigenvalues (1311) is marginally stable and its destabilization caused by dissipative, shear, 
and azimuthal magnetic field perturbations admits thus a natural interpre t ation as a 



dissipation-induced instability (IKrechetnikov fc Marsden 2007 



Kirillov 2007 



Kirillov 2009h . 



On the other hand, the matrix H can be considered as a result of a non- 



Hermitian complex perturbation Hi + H2 of a rea 
double semi-simple eigenvalues — diabolical points 



symmetric matrix, w hich has two 



Berry &: Dennis 20031 ) . This is a 



typical situation f or the problems of w a ve propagation in c hiral absorptive media 



(Keck et al. 2003: 


Berrv & Dennis 2003 


continua ( 


Kirillov 2009 


)• 



Kirillov et al. 20051 ) or in rotating symmetric 
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2.3. Dispersion relation in terms of dimensionless parameters 

The stability of the propagating plane wave perturbation fl25|) is determined by the 
roots 7 of the dispersion relation P(j) = det(H — jE) = 0, where 

P( 7 ) = 7 4 + a l7 3 + a 2 7 2 + (a 3 + *s)7 + a 4 + «&4 = 0. (32) 

We write the coefficients of the complex polynomial ( |32i) in the form 



a 2 = {oj v + lu v ) 2 + 2(oj a + lu„lu v ) + a 2 K 2 Q + Aa 2 uj\. 

a 3 = 2{u v + (jO v ){u a + lu v lu u ) + 2a 2 /to^ + 4a 2 (^ + uj u )uj a , 

a 4 = (a;^ + covu;,,) 2 — 4a 2 co>^fio + a 2 Kl(u\ + cj 2 ) + 4a 2 ^^cu 2 1 ^, 

6 3 = -8a 2 tt uj A uj A4t , 

b 4 = -Aa 2 n uj A uj A (2uj ri + uj u )-K 2 a 2 nQ l uj A uj A (uj ri -uj l ,). (33) 



After scaling the spectral parameter as 7 = X^/u! u uj v , we express the appropriately 
normalized coefficients (|33|) by means of the dimensionless Rossby number (Ro), magnetic 
Prandtl number (Pm), ratio of the Alfven frequencies (/?*), Hartmann (Ha*), and Reynolds 
(Re*) numbers 



_ 1 Ro dtt 
° ~ 20^ dR 



p v u v u A Q B° z 

, Pm = - = — , p = a — -, Re = a — , Ha = a- — . 
R=R rj u v u A lo v k^/nopvr] 



(34) 

Additional transformation yields the coefficients of the dispersion relation P(X) = in a 
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simplified form 

a i 

a 4 
h 
h 



2 ( v / Prri+ -^L= ) 



^ + 2(1 + Ha* 2 ) + 4/3* 2 Ha* 2 + 4Re* 2 Pm(l + Ro), 



ai(l + Ha* 2 ) + 2 ai /?* 2 Ha* 2 + 8Re* 2 (l + Ro)Vp^i, 

(l + Ha* 2 ) 2 + 4/T 2 Ha* 2 + 4Re* 2 + 4Re* 2 Ro(PmHa* 2 + 1), 



-8/3*Ha* 2 Re*v / Prri, 
-4/?*Ha* 2 Re*(2 + (1 - Pm)Ro). 



Therefore, we have exac tly reproduced the dispersio n relation of ( 



(35) 



Lakhin fc Velikhov 2007; 

— rr^ — r 

Riidiger fe Schultz 2008h . which generalizes that of flGoodman fc Ji 2002 ; iLiu et al. 20061 ). 



3. SMRI in the absence of the azimuthal magnetic field (j3* = 0) 



Let us first assume (3* = and study the onset of the standard magnetorotational 
instability (SMRI). The coefficients of the polynomial -P(A) are then real because 63 = 
and 64 = 0. We have 

1 



d\ = (Xi 



2 K/P 



m 



VP 



m 



"7 



a 2 = a 2 = j + 2(1 + Ha* 2 ) + 4Re* 2 Pm(l + Ro), 



a 3 = a 3 = a 1 (l + Ha* 2 ) + 8Re* 2 (l + Ro)v / Piri, 

a A = a 4 = (l + Ha* 2 ) 2 + 4Re* 2 (l + Ro(PmHa* 2 + l)). 



(36) 



Composing the Hurwitz matrix of th e real polynomial P(\) we write the Li enard and 



Chipart criterion of asymptotic stability (ILienard fc Chipart 1914 
A have ReA < if and only if 



Marden 19661 ): all roots 



d 4 > 0, d 2 > 0, hi — a% > 0, h 3 = aiaia 3 — d 2 d 4 — d 2 > 0. 



(37) 
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Explicit calculation of h% shows that it is a sum of squared quantities 

;2\2 /a2 



ha = 64 ^Pm*Re* 2 (Ro + 1) + ^\ + Ha* 2 a 2 + 4Re* 2 ^ > (38) 



Therefore, the condition /i 3 > is always fulfilled. 



The local definition of the Rossby number (1341) allows us to vary it for the background 
profile fl(R) = a + bR~ 2 changing the coefficients a and b because Ro = — b/{aR^ + b). On 
the other hand we can interpret the Rossby numbers as if they would correspond to quite 
general rotation profiles Q(R), which can have, e.g., the shape Q(R) ~ R w (with w = —3/2 
and Ro = w/2 = —3/4 for Kepler rotation). 

In the following we assume Ro > — 1 that corresponds to the centrifugally (Rayleigh) 
stable flow in the absence of the magnetic field. This reduces the conditions fl37|) to a 4 > 
that is equivalent to 



^l + Ha* 2 ) 2 + 4Re* 2 
4Re* 2 (PmHa* 2 + 1) 



Ro > R ° C = ~ L^^2 I TV ( 39 ) 



Note that in the absence of the magnetic field, Ha* = 0, the inequality ( |39l) is 

Ro > Ro v = -1 —J, (40) 

4Re* 

where we define the viscous Rayleigh line Ko v . In the inviscid limit Re* — > oo it is reduced 
to the Rayleigh's centrifugal stability criterion 

Ro > Ro 4 = -1, (41) 

where Ro* is the classical inviscid Rayleigh line. 

As is seen in the Figure [Deleft), there are two extrema of the function Ro c (Ha*) at 



Ha ,_ = Ml iWd-hJHP.y (42) 



which agrees with the results of (IJi et al. 200ll ). Triggered by the vertical magnetic field 



(B° z ^ 0) at some values of Ha* the flow becomes unstable for Ro > — 1. It can be stabilized 



-18- 



again, however, with the further in crease of Ha*, which is a hallmark of the standard MRI, 



cf. with Figure 1 in (IJi et al. 200ll ). 



The maximal values of the Rossby number at the peaks of the boundary of the SMRI 
domain are 

4Re* 2 Pm 2 + fPm - 1 + ^(1 - Pm) 2 + 4Re*' 2 Pm 2 Y 
j^ Q c V )__ (431 

4Re* 2 Pm 2 ^(1 - Pm) 2 + 4Re* 2 Pm 2 
The two extrema at Ha* max 7^ exist when the radicand in (|42|) is positive, that is when 

J (2 - Pm)Pm . , 

Re*p m > — J - . (44) 

Otherwise, the unique maximum is at the origin, Figure fright). This condition also follows 
from the positiveness of the second-order coefficient in the series expansion of Ro c (Ha*) at 
Ha* = 0: 

1+4R ^_2-(4R^ + l)P m 
4Re* 2 4Re* 2 V ; 

Setting Ro^ ax > —3/4, we find the conditions for existence of the standard MRI at and 
above the Kepler line in terms of the magnetic Reynolds number Rm* = Re*Pm 



Rm* > -Vl + 3Pm. (46) 
3 

For Pm <C 1 one should have Rm* > 2/3 to have SMRI for the Kepler flows, which leads to 
Re* ^> 1. At such values of Re* the following formal asymptotic expansions of Ha* max (Re*) 
and Ro max (Re*) are valid 

/ „ r 1 1 - Pm (1 - Pm) 2 , , 

Ha* max = ±V2Re* =f Ro c max = -— —— + = ? - \ 47 

4PmVRe* PmRe* 2Pm 2 Re* 2 2 3 Pm 3 Re* 3 V ' 

The asymptotic expansion (147p gives a simple scaling law which is known to be a 
characteristic of SMRI 

Re* = ^Ha* 2 . (48) 
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This equation is identical to N* := Ha* 2 /Re* = 2 where N* is often called interaction 
parameter (in technical magnetohydrodynamics) or Elsasser number (in geo- and 
astrophysics) . 

The SMRI can be interpreted as destabilization of slow Magneto-Coriolis waves 



flNornberg 2008 



Nornberg et al. 20091 ) . Indeed in the presence of shear, Ro 7^ 0, we find 
from the equation fl32|) with the coefficients fl33|) that in the absence of dissipation {uj v = 0, 
u v = 0) the eigenvalues are 



7 = ±y-2fijja 2 (l + Ro) -u 2 A ±2Q a^/Q 2 a 2 {l + Ro) 2 + u 2 A . (49) 
At the critical value Qq — where 

* = ~2a^MV < 50 > 

the branches of the slow magnetic-Coriolis waves merge with the origination of the double 
zero eigenvalue, see Figure [T^b). Splitting of this eigenvalue yields positive real eigenvalues, 
see Figure [Qc). Note that the SMRI threshold ( 1501) is equivalent to 

Ro = — 2 (51) 

4Re* 2 Pm v ' 

that follows from fl39|) when Ha* — > 00. 



4. HMRI in the presence of an azimuthal magnetic field {(3* 7^ 0) 



The fact that an additional azimuthal field ch anges the charac ter of the MRI drastically 
had been detected by Knobloch as early as 1992 (IKnobloch 19921 ). He had shown that in 



this case the instability appears in form of a travelling wave (see also ( Knobloch 19961 )). 
However, the difference in the scalin g behaviour for small Pm be tween standard and helical 



MRI was worked out only recently (IHollerbach fc Riidiger 20051 ). and is still the subject 
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of intense debate. In this section we will contribute to this discussion by focusing on the 
specific Pm dependence of the helical MRI. 



4.1. Bilharz criterion for asymptotic stability 



With the appearance of the azimuthal magnetic field {(3* ^ 0), the coefficients of 
the polynomial P(A) become complex. This breaks the symmetry of the eigenvalues with 
respect to the real axis of the complex plane and consequently may lead to dramatic changes 
in the stability properties of the system. 



In contrast to previous studies (ILakhin fc Velikhov 2007 



Riidiger k Schultz 2008T ) that 



were based on the study of approximations to the roots of the dispersion relation 

P(A) = A 4 + M 3 + a 2 A 2 + (a 3 + ib 3 )X + a 4 + ib 4 = 0, (52) 



we prefer to use the Bilharz criterion (jBilharz 1944 



Marden 19661 ) of asymptotic stability 



of the roots of complex polynomials. This criterion establishes the necessary and sufficient 
conditions for all the roots to be in the left part of the complex plane (ReA < 0) in terms 
of positiveness of the main even-ordered minors of the Bilharz matrix. For the polynomial 
P(A) with the coefficients (|35|) this matrix is 

/ , \ 



B 



a i 


-h 




















h 


a 3 




-h 














— a 2 





h 


a 3 


a,i 


-h 











— ai 


-a 2 





h 




a,\ 


-h 


1 








— a% 


-a 2 





h 










1 








-ax 


— a 2 

















1 








— a,\ 




















1 






(53) 
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The Bilharz stability conditions (IBilharz 1944 ; iMarden 19661 ) require positiveness of all 
diagonal even-ordered minors of B 

mi = a 3 a 4 + 6 3 o 4 > 0, (54) 

m 2 = (a 2 a 3 — a 1 a 4 )m 1 — 0364 > 0, 

m.3 = (oi«2 — 0,3)1712 — {<x{a±a2 + (01&3 — &4) 2 )mi + 0104(6402(264 — ai6 3 ) + 0^04) > 0, 

777.4 = ai?773 — aia 3 r?72 + (a| + a 2 04&3 — 2ai0 2 )mi + a^a^a^ — a 3 ) — o 2 a 2 a2 + h\ > 0. 



The inequalities (1541) determine the stability condition of the general dispersion relation 
(152]) in the presence of both vertical (B®) and azimuthal (-B9) components of the magnetic 
field. 

We first note that for /3* = the stability conditions (154]) are reduced to the stability 
condition a 4 > that was derived in the previous section. Indeed, with /?* = the 
coefficients 63 and 64 vanish to zero, which yields 

m! = 040.3, m 2 = 0.40.3(0-20.3 — 01(24), m 3 = 04(0-203 — 0104)^3, m 4 = a.4/13. (55) 

In view of S3 > and h 3 > it remains to check the sign of the expression 020-3 — 0104. 
Explicit calculation yields 

0203-0104 (4PmRe* 2 (Ro + l) + Pm + Ha* 2 + 2) 2 Pm 
2 PmyPm 

(4Re* 2 Pm 2 + l)(Pm + l)Ha* 2 +Pm 2 Ha* 2 (Ha* 2 + 3 + Pm) + 1 q 

PmvPm 

Therefore, for /3* = the conditions (155]) are reduced to the inequality a 4 > that 
determines the stability domain that is adjacent to the domain of SMRI. 

In Figure [3] we plot the boundary (1391) of SMRI domain to compare it with the domain 
of HMRI given by the inequalities (l54"j) . We see that the domain m 4 > is an intersection 
of all the domains > 0, i = 1,2,3,4. Thus, the region of HMRI, shown by dark grey 
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in Figure [31 is adjacent to the domain 1714 > 0. Although this fact is not a proof that the 
inequalities f[M|) are reduced to the last one, our numerical computations of the domains 
and of the roots of the dispersion relation as well as the analysis of the inductionless 
approximation in the next section confirm that = is the boundary of HMRI domain. 



4.2. Inductionless approximation 



As it was first observed in (IPriede et al. 2007I ). a remarkable feature of HMRI is that 



it leads to destabilization, even in the limit Pm — > 0, for some Ro > — 1, although not until 
the Kepler profile (Ro = —0.75). Below we prove this. 

Let us consider the Rossby number as a function of the magnetic Prandtl number and 
fix all other parameters. Substituting Ro = a + bPm + . . . into the equation m 4 = and 
collecting the terms with the identical powers of Pm, we find the quadratic equation on the 
coefficient a, which can be exactly solved. Therefore, in the limit Pm — > +0 there are two 
branches of the function Ro(Ha*, Re*, /?*): a positive one (Ro + > 0) and a negative one 
(Ro" < 0) 

Ro ± 

(1+Ha* 2 ) 2 +4/3* 2 Ha* 2 (l+/?* 2 Ha* 2 ) 



2/^Ha 



(57) 



± 



(2/5* 2 Ha* 2 +Ha* 2 +l)W(l+Ha* 2 )%4/3* 2 Ha* 2 (l+/3* 2 Ha* 2 )+2gi^ ((l+Ha* 2 ) 2 +4/3* 2 Ha*' 



2Ha* 4 /3* 2 

When (3* — > +0 the function Ro + tends to infinity while for Ro~ we get 

/3*^+o"~~ 4Re ' 



(1 + Ha* 2 )' 

lim Ro" = -1 - v — ^ 2 ; . (58) 



The expression (1581) can also be obtained as a limit of Ro c defined in fl39|) when Pm — > +0. 
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Calculating the derivative we find that it is strictly positive for all Re* G (0, +00) 

<9Ro~ 
<9Re* 

(1 + Ha* 2 ) 3 + /3* 2 Ha* 2 (6 + 8Ha* 2 + 2Ha* 4 + 8/3* 2 Ha* 2 ) 



2Re* 2 ^/Re* 2 ((l + Ha* 2 ) 2 + 4/3* 2 Ha* 2 (l + /5* 2 Ha* 2 )) + Ha* 4 /3* 2 ((1 + Ha* 2 ) 2 + 4/5* 2 Ha* 2 ) 
> 0. (59) 

Consequently, the maximal value of Ro~ < is attained when Re* — ► +00. In this 
limit the function Ro~(/3*) has a maximum 



r — (ri *\ y p — -Ha* 2 -2 + y/4 + 6Ha* 2 + 2Ha* 4 

ROmax(Ha ) = max hm Ro = ——2 (60) 

&* Re*-^+oo 2Ha* 

at 

ffi-CH.*) = V2 2 + H 2 " a ' 2 . (61) 
Since the derivative of the function Ro~ ax (Ha*) is strictly positive for all Ha* > 



dRo max (Ha*) 3Ha* 2 + 4-2V4 + 6Ha* 2 + 2Ha* 4 , . 
^ = 4 > (62) 

« Ha * Ha* 3 V4 + 6Ha* 2 + 2Ha* 4 

we conclude that the global maximum of the function Ro~(Ha*, Re*, j3*) coincides with the 
maximal value of Ro~ ax (Ha*), which is attained at Ha* — > +00 and is therefore 

max Ro" = maxRo max (Ha*) = 2 - 2^2 ~ -0.8284, (63) 

Ha*,Re*,/3* Ha* 

being exactly th e same value that was found in the highly resistive inviscid limit in 



( ]Liu et al. 20061 ). The corresponding optimal value of f3* in the limit Ha* — > +00 is 



/? max = ^ 0.7071. (64) 

Note that numerical maximization of Ro frequently leads to the extrema corresponding to 
the values of (3* ~ 0.7 even for Pm 7^ 0. 
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Extending the inviscid result of (ILiu et al. 20061 ) we establish that in the inductionless 
approximation (Pm = 0) the upper bound for HMRI is 



RcT(Ha*,Re*,/F) < 2 - 2^2 ~ -0.8284. 



(65) 



Proceeding similarly, we find that 



Ro + (Ha*, Re*, 0*) > 2 + V2 ~ 4.8284, 



(66) 



where the minimum is attained at the same extremal valu e of 0* given by (164"I) . The lower 
bound (166]) for Ro + exactly coincides with that found in (ILiu et al. 20061 ) in the highly 
resistive inviscid limit by analyzing the roots of the dispersion relation. However, it should 
be noticed that the character of this Ro + is still unclear. Since up to present we have not 
obtained any corresponding result from a ID global eigenvalue solver, it remains to be 
checked if this result is an artefact of the short wavelength approximation. 



Anyway, quite in accordance with (ILakhin fc Velikhov 2007 



Liu et al. 20071 ) we 



conclude that in the inductionless approximation there is no HMRI for 



2 - 2V2 < Ro(Ha*, Re*, 0*) < 2 + 2V2, 



which excludes HMRI for the Kepler law and for other shallower velocity profiles. 



(67) 



Finally, we would like to find a scaling law for HMRI to compare it with that of SMRI 
(j48|) . The HMRI scaling law for the maximum of the critical Rossby number at infinity 
(which works well, however, starting from Ha* ~ 3) reads 

V2 



Re* 



1 



Ha 



*3 



(68) 



In terms of the interaction parameter, this can be rewritten as N*Ha* = 1/(1 + ^/2). 
This scaling is rather different from the scaling of SMRI (Eq. 45). 
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4.3. HMRI in the case when Pm ^ 

In the previous section we have confirmed that for Pm = 0, HMRI does not work for 
Keplerian flows, at least according to the WKB approximation. Nevertheless, Hollerbach 
and Riidiger had shown that it does when considered as an eigenvalue problem, pr ovided 



that at least one of the radial boundary is conducting (IRiidiger fc Hollerbach 20071 ). 



In this section we analyze the dispersion relation without the simplifying assumption 
that Pm = 0. As it follows from the equation (1391) . in the absence of the azimuthal 
component of the magnetic field ((3* = 0) there is no SMRI for Pm < Pm c , where at the 
threshold 

(l + Ha* 2 ) 2 + 4Re* 2 (l + Ro) , . 

Pm = Pm c := — i '- = \ '-. 69 

4RoRe Ha 

The SMRI in this case develops when Pm > Pm c . In Figure [5] the threshold fl69|) is shown 
by the dashed line. 



4-3.1. First scenario of HMRI excitation 

When the azimuthal magnetic field is switched on ((3* 7^ 0), the instability threshold 
Pm c depends on (3*. At small values of (3* the threshold slightly increases so that 
the boundary of the domain of instability bends to the right of the dashed line in the 
(Pm, /3*)-plane, see Figure [51 The behavior of the instability boundary with further increase 
of (3* is determined by the Rossby number. 

For Rossby numbers close to Ro = — 1 the instability boundary at some (3* 7^ bends 
to the left, crosses the dashed line and forms a "semi-island" of instability with its center 
located close to f3* = Figure 0(a). This enlargement of the instability domain on 
the left of the dashed line is entirely a consequence of significant helicity of the magnetic 
field. Although the whole grey area of instability in Figure 0(a) is of course the region of 
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the helical magneto-rotational instability (HMRI), the main interesting effect of non-zero 
helicity is the excitation of instability at small and even infinitesimally small values of Pm 
in the range where SMRI is not possible. With the increase of Ro the specific effect of 
(3* becomes even more pronounced with the bifurcation of the instability domain to the 
isolated region (island) that lies on the left of the dashed line and to the "continent" on 
the right of it, Figure E(b,c). The island of HMRI illustrates both the destabilizing role of 
the azimuthal magnetic field component and the non-triviality of inducing HMRI below 
the threshold fl69l) at small negative values of the Rossby number (Ro > —1). For these 
reasons, we propose to refer to the instability for (3* ^ on the left of the dashed line as 
the essential HMRI and that on the right as the helically modified SMRI. 

Despite the apparent discontinuities in the (/?*, Ro)-plane, the three-dimensional 
domain of HMRI in the (/?*, Pm, Ro)-space has a smooth boundary given by the expression 
m-4 = 0, Figure[6^a). As it is seen in Figure[6](a), the function Ro c (Pm, f3*) has local extrema 
at some Pm ^ 0, yielding regions of HMRI that are separated from each other in the plane 
(Pm, /3*), Figure [5]^b,c). Since the maximum is attained at small but finite values of Pm, 
the corresponding boundary of HMRI in the (/?*, Ro)-plane at Pm ^ can exceed that in 
the inductionless limit (an instability induced by the viscosity u v ^ 0) and, moreover, the 
limiting bound Ro = 2 — 2v^2, as is clearly seen in Figure El The one-dimensional slices in 
the (/?*, Ro)-plane converge however to the region of the inductionless HMRI when Pm — ► 0. 
Therefore, in comparison to the inductionless limit, for Pm ^ we obtain higher values of 
the maximal Rossby numbers corresponding to the excitation of HMRI — a quite promising 
similarity of this local WKB a nalysis with the observation of HMRI for Keplerian flows 



with conducting boundaries in (IRiidiger fc Hollerbach 20071 ) . In Figure [7J the evolution of 



the stability boundaries in the (Pm, Re*)-plane with the increase of (3* demonstrates the 
details of the mechanism of reduction of the critical Reynolds number, which is another 
important characteristic of HMRI. 
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To clarify the nature of HMRI and SMRI and their relation to each other we inspect 
now the roots of the dispersion relation as functions of Pm. Series expansions of the roots 
in the vicinity of Pm = at j3* = yield 



Ai 



3 



A 



2.1 



-1 - Ha* 2 ± 2Re* v / -( 1 + Ro) + o(Pm 1/2 ), 

^ + Ha* 2 v / PnT + o(Pm 1/2 ) (70) 



VPin 

Therefore, two eigenvalues Ai^ branch from zero and the other two A 2 ,4 branch from minus 
infinity. The eigenvalues A 2 ,4 are real and negative, whereas the eigenvalues Ai^ are real in 
the vicinity of the origin for Ro < — 1 and complex otherwise with the frequency 

k 

io = 2ft -VRo + 1. 
k 

The eigenvalues Ai >3 correspond to inertial waves with uj = u c : = 2Q n k 7 Jk, if we assu me 



rotation without shear (Ro = 0) and without damping, see e.g. (INornberg et al. 20091 ). 



In the particular case, when (3* = and Ro = — 1 the dispersion equation is exactly 
solved 



-Pm - 1 ± \/{Pm - l) 2 - 4PmHa*(Ha* - 2VPni)Re* 
Al,2 = 



2VP 



m 



-Pm -1±\ (Pm - l) 2 - 4PmHa*(Ha* + 2 v / Pm)Re* 

A - = ~ • (71) 

The eigenvalues (17T|) are real near the origin, because Ro = — 1 



Ai, 3 = -v / P^-Ha* 2 v / P^±2Re*Ha*Pm + o(P 



m 



1 



A 2>4 = 7= + Ha* 2 VP^riT2Re*Ha*Pm + o(Pm). (72) 

V Pm 

Equating the first of the equations flTTj) or f!72|) to zero we reproduce the expression for the 
threshold (ED at Ro = -1. 
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In another particular case, when (3* = 0, Re* = 0, and Ro = the exact solution to the 
dispersion equation is two double eigenvalues 



A 



-- fv / Pni + —^J) ± \ - (VPm -Ha* : 



that are expressed in terms of the viscous, resistive, and Alfven frequencies in (130]) . 

Consider the eigenvalues corresponding to the values of parameters of Figure [5j Since 
Ro = —0.85 > —1, there are two real and two complex branches of eigenvalues A = A(Pm) 
at f3* — , see Figure [S](a). One of the real branches that comes from minus infinity 
changes its sign at the threshold ( l69i) and excites SMRI with zero eigenfrequency. Due to 
the large negative values of the real eigenvalues of the critical branch, there is no way to 
destabilize the flow at small Pm. However, new opportunities for destabilization occur with 
the increase of (3* that is accompanied by the qualitative change in the configuration of the 
eigenvalue branches. 

When the azimuthal component of the magnetic field is switched on, the real 
eigenvalues of the critical branch get complex increments, Figure El(b,c). With the increase 
of /?*, this critical real branch deforms and interacts with a stable complex one of an 
inertial wave until at (3* ~ 0.084 they merge at a point with the origination of the 



double comp 



f lBerrv 2004 



ex eigenvalue with th e Jordan block known as an exceptional point or EP 



Mailybaev et al. 20051 ). see Figure [8]^b). Notice that another exceptional point 
corresponds to negative j3* . With the further increase of (3* this configuration bifurcates 
into a new one, where parts of the stable and unstable branches are interchanged, Figure 
|8]^c). The new critical eigenvalue branch consists of c omplex eigenvalues that d emonstrate 



the typical generalized crossing scenario near an EP (jOr 1991 



Keck et al. 20031 ), when real 



parts avoid crossing while imaginary ones cross and vice versa Figure [S](c). 

In Figure [9] the " surgery" of eigenvalue branches is clearly seen in the complex 
(ImA, ReA)-plane. Although "on the surface" (for Re(A) > 0) nothing special happens, the 
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deep reason for the exchange of the fragments between the branches is "hidden" in the 
Re(A) < region at some finite value of Pm where an EP is formed, see Figure M^c). The 
critical branch that was responsible for SMRI leaves its stable "tail" coming from minus 
infinity (Figure [9](a,b)) and instead "catches" a fragment of a stable branch of complex 
eigenvalues with small real parts that comes from zero (Figure E2(d-f)). This re-arranged 
branch of complex eigenvalues is much more prone to instabilities at low Pm than the 
original critical one as the further increase of j3* confirms. Indeed, the negative real parts 
become smaller and around (3* = 0.7 there appears a new interval of HMRI at those values 
of Pm, at which SMRI did not exist, see Figure ITDT a). This interval exactly corresponds to 
the island of HMRI shown in Figure [5](b,c). We note here that the numerical calculation of 
the roots of the dispersion relation confirms the boundaries of the regions of HMRI given 
by the Bilharz criterion: m 4 = 0. 

The hidden exceptional point governs transfer of instability between the branch of 
(helically modified) SMRI and a complex branch of the inertial wave that after interaction 
becomes prone to destabilization. This qualitative effect explains why switching the 
azimuthal component of the magnetic field on we get HMRI as a travelling wave whereas 
SMRI was a stationary instability. Moreover, as Figure [TOT c) shows, the new critical 
branch is characterized by a broad band of unstable frequencies while the tail of the branch 
responsible for SMRI corresponds to a more sharply selected unstable frequency which is 
close to zero at (3* ^ 0. 

The above observations are in agreement with the observation of Liu et al. (2006) 
that, in contrast to SMRI, which is a destabilized slow Magneto- Coriolis wave, HMRI is a 
weakly destabilized inertial oscillation. Further results on interpretation of the HMRI as an 
unstable MHD-wave as well as on its relation to the dissipation-induced instabilities will be 



published elsewhere ( jFukumoto et al 



) 
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4-3.2. Second scenario of HMRI excitation 

The remarkable complexity of the phenomenon of HMRI manifests itself in different 
scenarios of destabilization. It turns out, that the transition from SMRI to HMRI through 
the exceptional point is not the only way to instabilities at low magnetic Prandtl numbers. 
At higher values of Re* and Ha*, in the presence of the azimuthal magnetic field the inertial 
wave can become unstable without mixing with the critical SMRI branch. 

In contrast to the scenario of the first type when one mixed complex branch becomes 
unstable at different intervals of Pm and causes both the essential HMRI and the helically 
modified SMRI, in the new situation the inertial wave causes the excitation of the essential 
HMRI and the critical real branch remains responsible for the helically modified SMRI, as 
is clearly seen in Figure [TTJ 

Most surprisingly, the inertial wave branch can become unstable twice with the increase 
of (3*. The first time this happens in the vicinity of (3{ = |\/2 ~ 0.707, see Figure HITb). 
then — in the neig hborhood of /3% = \y/2 ~ 2.121, as is visible in Figure [TTT c). In the 
(Pm, /?*)-plane this yields two islands of the essential HMRI that coexist with the continent 
of the helically modified SMRI, Figure fT2T a). The real parts of the unstable branches shown 
in black and grey in Figure [13] correspond to the first and second essential HMRI islands 
and to the continent of the helically modified SMRI, respectively. 

The difference between the two HMRI scenarios is visible also in the other parameter 
planes. For example, in (Ha*, Ro)-plane the domain of HMRI developed by the first scenario 
and corresponding to the HMRI-island of Figure fl2T c) has an SMRI-like form with two 
peaks, see Figure fT2T d). The second HMRI excitation scenario leads to the inclusions 
of the essential HMRI in the helically modified SMRI domain, Figure [T27 b). Although 
the islands in Figure fT27a) and Figure fT27c) look similar, the former are slices of the two 
three-dimensional domains that intersect each other along an edge, while the latter are 
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slices of the same smooth three-dimensional region of instability, cf. Figure H^a). 

5. Conclusions 

The helical magnetorotational instability is a more complicated phenomenon than the 
standard one. We found evidences that HMRI can be identified with the destabilization of 
an inertial wave in contrast to SMRI that is a destabilized slow Magneto- Coriolis wave. We 
established two scenarios of transition from SMRI to HMRI: the first one is accompanied by 
the origination of a spectral exceptional point and a transfer of instability between modes 
while in the second scenario two independent eigenvalue branches become unstable. We 
distinguish between the essential HMRI that is characterized by small magnetic Prandtl 
numbers at which SMRI is not possible, smaller growth rates than SMRI, and by non-zero 
frequencies and the helically modified SMRI which is caused by a small perturbation of 
the unstable real eigenvalue branch and is thus characterized by high growth rates, small 
frequency and relatively high magnetic Prandtl numbers within the usual range of SMRI. 
With the use of the Bilharz stability criterion we established explicit expressions for the 
stability boundary and proved rigorously the bounds on the critical Rossby number for 
HMRI in the inductionless limit (Pm = 0). Nevertheless, we revealed that for Pm ^ 
these bounds can be easily exceeded — an indicator in favor of the HMRI for small negative 
Rossby numbers. Finally, we found that for small negative Rossby numbers the essential 
HMRI forms separated islands that can coexist simultaneously in the (Pm, /?*) -plane. 

The work of O.N.K. has been supported by the research grant DFG HA 1060/43-1. 
The work of F.S. has been supported by the German Leibniz Gemeinschaft within its SAW 
programme. F.S. likes also to thank G. Gerbeth, R. Hollerbach, J. Priede and G. Riidiger 
for many stimulating discussions. 



-32 - 



A. Linearization with respect to non-axisymmetric perturbations 



We linearize equations ([I])-© in the vicinity of the stationary solution ([3J- (jTj) 
assuming general perturbations u = u + u', p = p + p', and B = B + B' and leaving 
only the terms of first order with res pect to the primed quantities. With the notation 



flGoodman fc Ji 2002 



Liu et al. 2006^ 



a d 
dt = d~r 



« 9 

OA, = — , 



9r= M , d z = -, 4 = ^+1 D = tf R d R +±- 2 dl+dl (Al) 
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we write the linearized equations in c ylindrical coordinates, cf. ( IGoodman &: Ji 2002 



Pessah fc Psaltis 2005 : 



Liu et al. 20061 ) 



(d t - vD)u' R - 2fitA + VLd.pu'n 



K 



(d t - vD)u'^ + — u' R + fld^ 



(d t - uD)u' g + ttd^ 



(d t - V D)B' R 
(dt-riD)B^ 



(d t - V D)B' z 





9 R P' + j Q {Bld R B' z + BpnB', - ^B'^j 



+ 



PoP 



B%B' R + ?±d,B' R - 2^5; 



u R 2 



1 1 



+ 



PoP 

1 



Po 

Bld^ + ^B'A v 



2 , ^ 



P 



+ 



Pop 
5° 



«V + - (B%B' Z + Bp z B'^) 
Po 

*M + §d,B' z 



-jjj-d+u'n - n^sj + B%u' R 



= B%u'^ - Bld z u' z - Bld R u' R - u' R d R B° 4 
+ uld R B' R + u°d z B' z + B' R d R u° 



B ( : 



B° 



-B° z d R u' R -^u' R -^- 

9 R U R + ft 9 <l> U '<f> + 



^B' z + ^u' z 
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a) b) c) 



Fig. 1. — (a) Coriolis splitting of the Alfven plane wave into the fast and the slow Magneto- 
Coriolis (MC) waves, (b) shear causes interaction of slow MC branches with the origination 
of the double zero eigenvalue, (c) splitting of the double zero eigenvalue yields positive real 
eigenvalues (SMRI). 
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Fig. 2. — The domain according to equation (J39l) of the standard MRI (grey) for Pm = 10 5 
and (left) Re* = 2 • 10 5 or (right) Re* = 0.5. The left plot shows the typical SMRI peaks 
flJi et al. 20011 ). 
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Fig. 3. — The domain according to equation ( l39|) of SMRI (light grey, (3* = 0) and that of 
HMRI (dark grey) for Pm = 10~ 5 and Re* = 2 • 10 5 . The thin full line marks the stability 
boundary Ro = Ro c . The dashed, dotted, and bold lines bound the domains where the 
Bilharz minors mj .2,3,4 > 0. The domain of HMRI which is adjacent to the intersection of 
these domains is shown in dark grey. The boundary to the HMRI domain (bold line) is 
m 4 = 0. 
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Fig. 4. — Inductionless approximation (Pm = 0): Two domains (1541) of HMRI (grey) are 
within the bounds (1671) . It is seen how the scaling law fl68|) works. HM RI is possible even for 
Pm = - the paradox of inductionless magnetorotational instability (IPriede et al. 20071 ). 
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a) b) c) 

Fig. 5. — The relation of SMRI (right to the dashed line for f3* = 0) and HMRI (grey area). 
For (3* a (semi)-island of HMRI exists in the neighborhood of (3* = y/2/2 for the values of 
Pm at which SMRI was not possible in the absence of the azimuthal magnetic field (/3* = 0). 
The size of the HMRI island decreases with the increase of Ro. 




Fig. 6. — (a) The three-dimensional HMRI domain shows how the slices of it in the plane 
Pm = const converge smoothly to the inductionless HMRI domain for (3* = 0. (b,c) HMRI 
domain (7B4 > 0) can exceed the inductionless bound Ro = 2 — 2^/2 when Pm 7^ 0. Thin lines 
show the boundary of the inductionless HMRI domain, bold ones show that with Pm 7^ 0. 
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Fig. 7. — SMRI boundaries and those of HMRI in the plane (Pm, Re*) for various val- 
ues of Ro. A considerable deformation of the boundaries when (3* is switched on is 
seen in the region of sm all Pm. The pictures are in agreement with the calculations of 
jRiidiger fc Schultz 20081 ). 
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Fig. 8. — Ha* = 5, Re* = 100, Ro = —0.85: Real and imaginary parts of eigenvalues A 
as functions of Pm. (a) Two stable complex branches (inertial modes) and two real ones, 
one of the latter (grey) becomes positive and causes SMRI; (b) merging of the inertial mode 
with the deformed unstable SMRI-branch (grey) with the origination of an exceptional point 
(EP); (c) bifurcation yields new mixed eigenvalue branches that create new opportunities 
for instability at low Pm. 
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Fig. 9. — Ha* = 5, Re* = 100, Ro = —0.85: Eigenvalue trajectories in the complex plane 
when Pm changes from zero to Pm = 0.011 show the effect of transfer of instability between 
branches through a spectral exceptional point (EP). (a) SMRI is caused by the transition of 
a pure real eigenvalue from negative values to positive, (b) deformation of the critical branch 
(grey) for (3* ^ 0; (c) merging of the critical and stable branches with the origination of the 
double complex eigenvalue with the Jordan block - exceptional point (EP); (d-f) splitting 
of the EP and bifurcation of the eigenvalue trajectories: after a surgery the stable complex 
branch acquires an unstable tail of the critical branch (grey) and becomes unstable near the 
origin, which corresponds to an unstable traveling wave of HMRI. 
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Fig. 10.— Ha* = 5, Re* = 100, Ro = -0.85, (5* = 0.7: (a) Real part of the critical eigenvalue 
branch is positive while Pm is within the island or the continent of instability; (b) imaginary 
parts of stable and unstable branches; (c) trajectory of the critical eigenvalue that creates 
instability with the change of Pm from zero to Pm = 0.011. 
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Fig. 11.— Ha* = 17, Re* = 28850, Ro = -0.88: (a) Stable inertial waves (black) and 
critical SMRI branch of real eigenvalues (grey) for (3* = 0; (b) for f3* = 0.7 one of the inertial 
wave branches becomes unstable independently of the SMRI branch; (c) for (3* = 2 the same 
inertial wave branch becomes unstable a second time at higher frequencies independently of 
the SMRI branch. 
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Fig. 12. — (a) The first and second islands of the essential HMRI developed according to 
the second scenario in the neig hborhood of p* = ^2/2 and (3* = 3^2/2, (b) an essential 
HMRI inclusion in the helically modified SMRI domain, (c) an island of the essential HMRI 
developed according to the first scenario, (d) the corresponding view in the (Ha*, Ro)-plane. 
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Fig. 13.— Ha* = 17, Re* = 28850, Ro = -0.88: (left, (3* = 0.7) Real part of the first critical 
eigenvalue branch (grey) is positive when Pm is within the continent while for the second 
critical branch (black) it is within the first HMRI island; (right, (3* = 2) real part of the first 
critical eigenvalue branch (grey) is positive when Pm is within the continent while for the 
second critical branch (black) it is within the second HMRI island. 



